Quantum Noise as an Entanglement Meter 
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Entanglement entropy, which is a measure of quantum correlations between separate parts of a 
many-body system, has emerged recently as a fundamental quantity in broad areas of theoretical 
physics, from cosmology and field theory to condensed matter theory and quantum information. 
The universal appeal of the entanglement entropy concept is related, in part, to the fact that it 
is defined solely in terms of the many-body density matrix of the system, with no relation to any 
particular observables. However, for the same reason, it has not been clear how to access this 
quantity experimentally. Here we derive a universal relation between entanglement entropy and 
the fluctuations of current flowing through a quantum point contact (QPC) which opens a way 
to perform a direct measurement of entanglement entropy. In particular, by utilizing space-time 
duality of Id systems, we relate electric noise generated by opening and closing the QPC periodically 
in time with the seminal S = h log L prediction of conformal field theory. 



Recent years have witnessed a burst of interest in the 
phenomena of quantum entanglement, and in particular, 
in entanglement entropy, a fundamental characteristic 
describing quantum many-body correlations between two 
parts of a quantum system. This quantity first emerged 
in field theory and cosmology [l|, Q under the name 
of "geometric entropy," and subsequently was adopted 
by quantum information theory. The notion of entan- 
glement entropy has provided a framework for analyz- 
ing quantum critical phenomena [3, 0, [H| and quantum 
quenches @, Q B 0- Recently it was used as a probe 
of complexity of topologically ordered states [lOj, [ill, [ij] ■ 
In addition, this quantity is of fundamental interest for 
quantum information theory as a measure of the re- 
sources available for quantum computation [l^] as well 
as for numerical approaches to strongly correlated sys- 
tems 



Can the entanglement entropy be measured? Here we 
identify a system where such a measurement is possible, 
thereby offering an affirmative answer to this question. 
In particular, we establish a relation between the entropy 
and quantum noise in a quantum point contact (QPC) 
[isl , an electron beam-splitter with tunable transmission 
and reflection. In essense, the QPC serves as a door be- 
tween electron reservoirs, which can be opened and closed 
on demand (see Fig[T]) . We show that the fluctuations of 
electric current flowing through the QPC can be used to 
quantify the entanglement generated by the connection, 
and thereby measure the entanglement entropy. 

On some level, the very idea of measuring a quantity 
that encodes information about many-body correlations 
of a large number of particles, which is what the entan- 
glement entropy is, may seem totally bizarre. Yet, as we 
shall see, the situation with the entanglement entropy is 
different from, for example, the many-body density ma- 
trix that depends on coordinates of all particles in the 
system and is thus indeed very difficult to measure. In 
the free fermion QPC problem analyzed below, all multi- 
particle correlations in the Fermi sea that are relevant for 
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FIG. 1: Schematic of a quantum point contact (QPC) with 
transmission changing in time. The left and right leads are 
initially disconnected, then connected at to < t < ^i, and 
then disconnected again. Electron transport, taking place at 
to < t < tl makes electrons delocalized among the two leads, 
generating entanglement and current fluctuations. 



entanglement are fully accounted for by temporal corre- 
lations of current ffowing through the QPC. As a result, 
perhaps somewhat surprisingly, noise measurement pro- 
vides sufficient information needed to determine the en- 
tanglement entropy. 

The situation considered here, which involves connect- 
ing and then disconnecting two parts of the system, is 
in a sense dual to the more conventional picture used in 
Refs.[3,[i,i,[l. In the latter approach, the many-body 
ground state of a translationally invariant system is ana- 
lyzed using a finite region in space. In our case, a window 
in time is used, to < ^ < ^i, during which particles can 
delocalize among the reservoirs, making them entangled. 

The relation between en tang lement and noise is at the 
heart of recent proposals to use current parti- 

tioning by scattering on a QPC for producing entangled 
particle pairs. Ref. [l9| considered measuring the entan- 
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glement entropy, emphasizing its relation with particle 
number statistics, however without making connection 
to electron transport and current fluctuations. We also 
note that generation of entanglement was recently ana- 
lyzed for critical Hamiltonians , for generic Hamiltoni- 
ans 0,01, as well as in a QPC under bias voltage [20| . 

Our analysis reveals a relation between entanglement 
production in a driven QPC with the theory of quan- 
tum noise [3], describing quantum and classical current 
fluctuations in this system. This approach, known as 
Full Counting Statistics (FCS), describes the probabihty 
distribution of transmitted charge using the generating 
fimction x(A) = X]^-cx3 Pne^'^^, where P„ is the proba- 
bility to transmit n charges in total. The function x(A) 
encodes all FCS higher moments, or cumulants Cm- 



logx(A) = J2 

m—1 



(1) 



where Ci, C2, C3... describe properties of the distribu- 
tion Pn such as the mean n, the variance {{n — fi)^), the 
skewness {{n — n)^) , etc. The cumulant C2 is available 
from routine noise measurement. Recently, C3 has been 
measured in tunnel junctions [U, and in QPC [2^, 
while cumulants up to 5th order where measured in quan- 
tum dots jM, i25,] . 

Below we establish a universal relation between FCS 
([1]) and the entanglement entropy generated in the pro- 
cess shown in FigU] Only even cumulants are shown to 
contribute to entropy: 
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Be = j2---)- The first few contributions are: 
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(3) 



This relation, which is completely general and valid for 
arbitrary driving, can be used to determine the entan- 
glement entropy from measured values of FCS moments. 
Similar relation can be derived for other quantities of in- 
terest, such as Renyi entropies, or single copy entropy 
23. 

In particular, quantum noise generated in the QPC 
switching on and off (Fig[T]) directly corresponds to the 
entanglement entropy found for conformal field theory, 
S — ^ ^ogL, where L is the size of window in space and 
c the conformal charge 0, H, . In this case the current 
fluctuations are gaussian {Cm^2 ~ 0), with a logarithmic 
variance C2 = \ log , where r is a short time cut- 
off set by the QPC switching rapidity. Combined with 
Eq.(l3]) this gives entropy S ^ \ \o'g\t\ — to|- Below we 
discuss how this logarithmic dependence can be verified 
using the setup shown in Fig[TJ 



Entanglement entropy is conventionally defined as the 
von Neumann entropy Si^p) = — Trplogp, where p is the 
reduced density matrix of a pure quantum state, made 
"impure" by confining it to a certain space region [l|, 
01 . In our case, the many-body state evolves as a pure 
state while the QPC is open (see FigH]), after which the 
reduced density matrix of the lead L is given by 



PL(ii) = T>;^(U(ii,toVoUt(to,ii)). 



(4) 



Here po is the initial density matrix of the system, U is 
the many-body evolution between tg and ii, and Tr^ is 
a partial trace over degrees of freedom in the lead R. 

Entropy production in the lead L as a result of QPC 
opening and closing is described by 



A5 = 5(pL(i))-5(po) 



(5) 



where the last term accounts for the entropy in the initial 
state. Because at finite temperature both terms in Eq. ^ 
are proportional to the lead volume, they can be large for 
macroscopic leads. The increment A5, however, remains 
well defined regardless of the lead volume. 

Below we focus on the zero temperature case when po 
is a pure state, described as a filled Fermi sea in the full 
system L + R, in which case the second term in Eq.© 
vanishes, giving AS = S{pL{t)). We associate with po a 
Fermi projection operator n in the single-particle space 
{E\n\E') = SE.E'diEp - E), where Ep is the Fermi en- 
ergy. The evolved system is described by a rotated Fermi 
projection n(j — UnU\ where U is the unitary evolution 
of the single-particle modes. 

Our first step will be to express the entropy in terms 
of single-particle quantities. For a generic gaussian state, 
Wick's theorem for operator products is satisfied in L+R, 
and therefore in particular it holds in L [sij. Therefore 
the reduced density matrix p^ is also gaussian: 



PL 



in L) 



(6) 



for some H . We define a single-particle quantity rriij — 
TvpLa\aj. For the evolved system, described by n(7, 
Wick's theorem gives rriij — {nij)ij- In what follows it 
will be convenient to extend m to L + R by setting 

M^PluuPl, (7) 

with Pl a projection on the modes in L, so that M = m 
in L and AI — in R. 

Entropy can be expressed through rriij for a generic 
gaussian state ©. Because of Fermi-Dirac statistics, 
TO = (1 + e^)^^, which gives H = log(m~^ — 1) [sij. 
Extending to to M in L + R, we write the entropy as 

S{p) = -Tr [M log M + (1 - M) log(l - M)] (8) 

where now the trace is taken in the space of single-particle 
modes in L. 
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Transport in a QPC is described by time-dependent 
transmission and reflection amplitudes A{t), B(t). In a 
Schrodinger representation, the scattering states are 
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\x) 







X <{) < x{t), 



and U\x)l,r = \x)l,r otherwise. Here x{t) — x + vpt, 
tx — —x/vp is the time of arrival at the scatterer, w^? is 
the Fermi velocity, and \x) l.r describes incoming (x < 0) 
and outgoing {x > 0) wavepacket states in the leads. 

In FCS approach it is convenient to work in a time 
representation [s^l , labeling states by times of arrival at 
the scatterer tx- In this representation the initial Fermi 
projection is given by n{t,t') = 2^i(t-t' +io) ^ ^^^^ ^ ^ 
2x2 identity matrix (in L, R). The evolved state njj is 



nu{t,t') = U{t)n{t,t')U\t'), U{t) = 



B{t) A{t) 
-A{t) B{t) 



These relations are inessential for the derivation of our 
main result, they serve here illustration purpose only. 

Our next step is to relate the quantity M, Eq.d?]), and 
the FCS generating function ([1]) which can be expressed 
as a functional determinant [16]: 



X(A) = dct(l -n + 



(9) 



This determinant must be prop erly regularized for in- 
finitely deep Fermi sea [H, [2g|. For our purposes we 
proceed to treat it as a finite matrix although a more 
rigorous treatment using C* algebra techniques is possi- 
ble In the spirit of 0,[2j, we rewrite jS]) as [30]: 



det(e*^-P^"(l -n + 



-i\PL 



Using the identities e 



(1 — n) = 1 — n and 



giAPi,(i n)^ — n we find: 

X(A) = det((l - n + nWe'^P^U)e-'^P^'^). (11) 

Next, we insert U'^U in the determinant to obtain 

X(A) = det(;7(l 

= det((l -nu + nue'^P^)e-'^^^^''^<^) 
= det((l + nuPLie'^ ~ l))e-*^(-P^")^), 



nWe'^P^U)e-'^P^''W) (12) 



where in the last line we used: e' 
Finally, noting that det(l + APl) 
any matrix A, we arrive at 



= 1 + (e*^ - 1)Pl. 
dct(l +Pl APl) for 



X(A) = det ((1 - M + Afe*^)e-*^("^^)^ 



(13) 



where M is the quantity ([7]) which defines the entropy. 

Now, with the help of the relation we can ex- 
press the spectral density of M, Eq.([7]), which lies be- 
tween and 1, through x(A). Indeed, changing param- 
eter A to z = (1 
minant: xi^) = 



yields z 

z 

From that we can write the spectral density of M as 
1 



M under the deter- 
det ((z - M)e-'("^^)^^W(1 - e*^(^))). 



where the coefficients A and B depend on dimM, 
Tr(nPi)[/, and Ci. Hereafter we ignore the delta func- 
tion terms because z — 0, 1 do not contribute to the 
expression ([8]) for the entropy which we rewrite as 

S = - [ dzn{z) [z log z -I- (1 - z) log(l - z)] (15) 
Jo 

Now it is straightforward to evaluate S by substitut- 
ing (P) into (Il4l),(IT5l) with A(z) = tt - i log (i - l) , 
and integrating by parts over z. We obtain series S = 
J2m=i ^^"i' where the coefficients am, after changing 
variable to u = i \og{j^), take the form 



(-2)^ 



1 U ITT 

, cosh u V 2 



(16) 



Finally, for m > 2, shifting the contour of integration as 
u u — i^, yields [32,] our main results ^ and ([3]). 

In order to compute FCS, it is convenient to return 
to the expression Q and apply the Riemann-Hilbert 
(RH) method introduced in [2§]. In this approach, 
one must factor the time-dependent matrix R{t) 



iXPL 



in ([9]), i.e. find matrix valued func- 



tions X±{z), analytic in the upper/lower half plane of 
complex z, respectively, such that on the real line: 



(17) 



with normalization X±{z) — > / at \z\ oo. 

We consider QPC switching between the on and off 
states several times < t[^^ < ... < < t[^K The 
RH problem is solvable in the case of abrupt switching 
because R commutes with itself at different times: R — I 
in the off state, R = e-^^^^s^^^-) in the on state. The 
solution of the RH problem is then given by the functions 



fx, z-4*'±i0' 



,(18) 



To find the determinant ^ with these X± we use the 
RH method [l^ to evaluate the derivative of logx: 

logx(A) = y tr (^±,X^^dtX^dxR^ dt - -^G, 

(19) 



G 



N 



log 



A^) _ AJ) Ai) _ Aj) 



Ai) Aj) ' 



where for i = j the denominators must be replaced by a 
short-time cutoff r. This gives gaussian charge statistics 



X(A) 



C2 - TT^G. 



(20) 



Because in this case the only nonvanishing cumulant is 
C2, we have S = ^C2, which gives S = ^ log for a 
^(z) = ^Im^^ logx(2; - iO) + AS{z) + B6{z - 1), (14) gj^gle QPC switching, in agreement with the 5 - i logi 
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Pulse width w/T 

FIG. 2: Entanglement entropy and current fluctuations. 
Noise power in a QPC driven by a pulse train, plotted 
as a function of the pulse width. Parameters used: driving 
frequency v = 500 MHz, short-time cutoff r = 20 ps. The 
noise as well as the entropy production are symmetric under 
w — > T — w. Note that at a narrow pulsewidth w <^ T the 
dependence reproduces the | log L behavior of the 
entropy. 
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FIG. 3: The effect of imperfect transmission in QPC. Spectral 
measure (|24|) scaled by G/27r^, for a QPC driven by pulses 
with transmission D in the "on" state (see Fig[21). At D < 1 
the entropy p5[) is reduced by a constant factor (inset) with 
the S logsiuTTi^w dependence unchanged. 



relation The case of multiple switching provides re- 
alization of the situation studied in Ref. . 

These predictions can be tested by measuring noise in 
a QPC driven by a periodic train of pulses (see Fig. 2). 
For N identical pulses, the relation at large N yields 



C2{N) 



N sin TTvw 



V = l/T. 



(21) 



Thus in a periodically driven QPC there is a fi- 
nite entropy production per cycle at a rate dS/dt = 
ii^log2i|f^. Fluctuations ^ with Ca cx TV corre- 
spond to electric noise with spectral power 



So 



log- 



(22) 



dependence ([^ becomes S2 = ^ log ^ , identical to 
the entropy for a single pulse. 

The result ((22)) must be compared with thermal 
noise. At a driving frequency v — 500 MHz, the ef- 
fective temperature corresponding to (I22p is Tcs — 
_hj^ log « 25 mK. In practice it may be possible 

to relax the constraint due to small Toff by detecting the 
noise at frequencies somewhat higher then ksTcs/h, 
detuned from the thermal noise spectral window. 

How sensitive are these results to imperfections in QPC 
transmission? It is straightforward to incorporate trans- 
mission D — \ A\^ < 1 in the "on" state in the RH anal- 
ysis because the matrices R{t) in ([9]) still commute at 
different times. Instead of e^'"^, the eigenvalues of the 
R matrix are now e**^* 
Making this change, we obtain 



with sin^A, = ^/Dsin^X [32 



X(A) 



tG 



(23) 



with G given by as above. Because this xW is non- 
gaussian, with nonzero higher cumulants, the simplest 
way to find the entropy is to use its relation with the 
spectral density of M, Eqs.^,^. Using ^ along 
with the relations between A,, A and z, we find 



G D 



Re 



Il-2z| 



7r2z(l-z) y/D^ -ADz{l- 



(24) 



at frequencies below v. For a short pulsewidth w, the 



As illustrated in Fig. 3, at 13 < 1 the function /i(z) van- 
ishes in the interval z„ < z < z+, z± = i (l ± \/l — . 

The entropy, found from (fT5|) and ([24]) . has the same 
logarithmic dependence p9|) on the times tg \ as above, 
albeit with a D-dependent prefactor. Thus, the predicted 
dependence S ~ logsinTrz/w remains robust. The behav- 
ior of the rescaling factor F = S{D)/S{1) (Fig. 3 inset) 
indicates that entropy reduction due to imperfect trans- 
mission in QPC can be attributed mostly to the change 
in the second cumuiant, G2 = -^G, with a relatively 
small correction due to higher cumulants. 

From the quantum information perspective, it is im- 
portant to isolate the part of the entropy accessible to 
local operations (i.e. respecting particle conservation in 
each lead) [23|. The particle-number restricted entropy 
may be easily obtained from xW- For the gaussian case 
considered above the change in the entropy due to this 
restriction is inessential (to be published). 

In summary, we have shown that the entanglement en- 
tropy can be directly inferred from statistics of current 
fluctuations. We derived a general relation between the 
entanglement entropy and electron transport via the full 
counting statitics. Builiding on this universal relation, 
we propose noise measurement in a QPC as a way to 
test theoretical predictions for the many-body entangle- 
ment in a realistic setting. This provides a new method 
to invetstigate many-body entanglement, and in particu- 
lar, its generation in non-equilibrium quantum systems. 
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